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RESPIRATION AND HEMOLYSIS 
OF ULTRAVIOLET 
IRRADIATED FROG ERYTHROCYTES’ 


Lewis E. Harvie AND S. P. MARONEY, JR. 
University of Virginia 
INTRODUCTION 


Ionizing radiations have been known for some time to cause hemo- 
lysis in mammalian erythrocytes. Early experiments indicated that X- 
irradiation-induced hemolysis of human and beef erythrocytes was 
preceded by slow swelling of the cells (Ting and Zirkle, 40). In the 
same report, these authors suggested that hemolysis and swelling were 
santa on an increased permeability to sodium and potassium ions. 


Subsequent analysis by Sheppard and Beyl (’51) confirmed the hypo- 
thesis that the cell was more permeable to sodium and potassium after 
X-irradiation; they further showed that these ions exchange quantita- 
tively. Non-ionizing radiations have also been shown to cause hemolysis 
(Sonne, ‘29; Cook, ’56; Green, ’56). The hemolytic effects of ultraviolet 
light as reported by Cook (56) and Green ('56), are, in general, in 
agreement with those produced by X-rays. 


It has been demonstrated that the energy for maintaining the 
balance of potassium and sodium in non-nucleate erythrocytes is supplied 
by glycolysis (Maizels, °51; Weller and Taylor, ’51; Solomon, °52). 
Green (56) has shown, however, that ultraviolet irradiation of human 
red cells causes an acceleration in the loss of potassium and gain of 
sodium ions without inhibiting glycolysis, Sheppard and Stewart (’52) 
obtained similar results with human red cells exposed to X-rays. This 
suggests that, in addition to metabolic energy, other factors are impor- 
tant in regulating the exchange of these cations in non-nucleated red 
cells. One obvious structure for such regulation is the cell membrane 
and, indeed, gross changes in the mammalian erythrocyte membrane 


This paper was presented to the Virginia Academy of Science on May 8, 1959 at the 
University of Virginia. An abstract has been published in the Virginia Journal of 
Science, Sept. 1959, 10(4): 263-264. 
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have been observed following X-irradiation (Zacek and Rosenberg, 50). 
It seems clear, then, that radiation-induced hemolysis in the mammalian 
erythrocyte is preceded by a change in ion permeability which is not 
dependent on an alteration in the energy supply of the cell, ie. 
glycolysis. 


Most of the investigations on erythrocytes have been concerned with 
mammalian cells; comparatively few studies on nucleated erythrocytes 
are available. Hunter and Hunter (57) have shown that the oxygen 
consumption of nucleated erythrocytes is significantly higher than that 
of non-nucleated erythrocytes when compared at the same temperature. 
Evidence for an active, oxidative metabolism in nucleated red cells of 
the chicken may be found in the work of Rubinstein and Denstedt (’53), 
and Defendi and Pearson (’55). Ion transport in nucleated erythrocytes 
(chicken, grass-snake) depends on this aerobic respiration (Maizels, 
54). The problem, then, arises as to the effect of hemolytic doses of 
radiation on the aerobic metabolism of the nucleated red cell. 


Martin and her co-workers (’56), in studying bullfrog erythrocytes, 
observed an inhibition in oxygen consumption of these cells following 
X-irradiation, The present work has been undertaken to study the rela- 
tion between ultraviolet-induced hemolysis and respiration following ul- 
traviolet radiation in nucleate erythrocytes of the leopard frog, Rana 
pipiens. The results indicate that this radiation increases hemolysis with- 
out inhibiting respiration. 


METHODS 


Frogs used in this investigation were kept in covered porcelain 
ans in about an inch of water, the temperature of which varied between 
20°C and 22°C. 


Blood from the aorta of one or two pithed frogs was drawn into a 
fine glass pipette, the tip of which had an inside diameter of 0.7 mm, 
and mixed with phosphate buffered Ringer’s (0.11M NaCl, 0.02M KCl, 
0.01M CaCl,, 0.0009M Na:HPO,, and NaH,PO,; pH 7.4) in which 
a few crystals of sodium citrate had been dissolved as an anticoagulant. 
Just before drawing up the blood, the pipette was rinsed with some of 
this solution in order to prevent clotting of the blood on contact with 
the inside wall of the pipette. Erythrocytes harvested from this blood by 
centrifugation were washed twice and finally suspended in phosphate 
buffered Ringer's solution, Cell concentration, determined from hemacyto- 
meter counts, ranged from 70,000 to 150,000 per cubic mm. 


The cell suspension was divided into two 3.5 ml aliquots, each of 
which was placed in a 50 ml round bottom quartz flask and rotated at 
five rpm to keep the cells in suspension. One aliquot was irradiated 
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with a general Electric G15T8 germicidal lamp, the times of irradiation 
ranging from 10 minutes to 60 minutes. The distance of the surface of 
the cell suspension from the lamp was 4.5 cm. All experiments were 
carried out at room temperature (approximately 25°C). The second 
aliquot, which was treated in the same way as the first except that it 
received no ultraviolet radiation, acted as the control. Respiration of 
both irradiated and non-irradiated cells was measured using a standard 
Warburg technique, following the procedures described by Umbreit 
et al (’57). Flasks containing two ml of cell suspension, to which was 
added glucose in a final concentration of 50 mg %, were placed in a 
water bath at 27.3°C and oscillated 100 times per minute through a five 
em stroke. Since cell number varied with each experiment, oxygen 
consumption was recalculated on the basis of 100,000 cells. 


Separate flasks containing about one ml of cells were oscillated in 
the water bath. For hemolysis determination, samples were taken from 
these flasks, diluted 1:10 with Ringer’s, and placed on glass slides for 
counts of hemolysed cells. Percent hemolysis was calculated as the 


number of hemolysed cells per total cells counted (total cell counts 
between 400 and 500). 


RESULTS 


Non-hemolysed cells had an opaque, yellowish appearance, with 
the nucleus clearly discernable. Dark bands running at right angles to 
the Iong axis of the cell were apparent. These have been described by 
Trotter (56) who attributed them to the residual alkali in the glass 
slide. As the cell swelled following irradiation, the bands disappeared, 
the cell became subspherical and at hemolysis rapidly faded with the loss 


of hemoglobin to a “ghost.” The appearance of the “ghost” was the end 
point in these experiments. 


Cumulative oxygen consumption increased as a direct linear function 
of time for about the first six hours from the beginning of the mano- 
metric measurements. A typical experiment (30 minutes ultraviolet 
radiation) is presented in Figure 1. It should be noted that there was 


no significant difference in the oxygen consumption of irradiated and 
non-irradiated cells. 


Hemolysis of irradiated erythrocytes increased in a similar manner 
for about the first four hours, This was followed by a leveling of the 
hemolysis curve. This “plateau” in the curve lasted three or four hours 
after which it rose sharply to 100%. Time required for complete hemo- 


lysis was approximately twelve hours. Non-irradiated cells showed no 
hemolysis. 


The initial amount of hemolysis after irradiation varied according 
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FIGURE 1 
Oxygen consumption and percent hemolysis of ultraviolet irradiated and 
non-irradiated frog red cells (open circles = respiration, ultraviolet; 
closed circles = respiration, control; open triangles = hemolysis, ultra- 
violet; closed triangles — hemolysis, control.) Ultraviolet dose = 30 
minutes. 
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PERCENT HEMOLYSIS 


FIGURE 2 


Respiration rates and per cent hemolysis values of ultraviolet irradiated 
and non-irradiated frog red cells four hours after exposure to radiation 
(open circles = _ respiration, ultraviolet; closed circles — respiration, 
control; triangles = hemolysis). The control respiration rate is included 
for each ultraviolet dose to better point up the fact that there is no 
significant difference in respiration rate between irradiated and non- 
irradiated cells. 
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to the dose used. At a 10 minute dose, no hemolysis was observed in 
the irradiated cells until three or four hours following irradiation. At a 
15 minute dose, hemolysis began about one hour after irradiation, 
while for a 20 minute dose, hemolysis began immediately following irra- 
diation. At 30 and 60 minute doses, 4 and 10 percent, respectively, of 
the cells were hemolysed in the first count following irradiation. 


Respiration rate for the first four hours and the per cent hemolysis 
at four hours following irradiation are presented as functions of ultra- 
violet dose in Figure 2. The per cent hemolysis increased with the 
ultraviolet dose while the respiration rate was not significantly different 
between the irridiated and non-irradiated cells. 


DISCUSSION 


The first apparent change following ultraviolet radiation was a 
swelling of the frog red cell. A reasonable explanation for this swelling 
involves the phenomenon of colloid osmotic pressure which consists es- 
sentially of two phases (Ponder, “48): (1) a lytic agent alters the 
cell membrane’s selective permeability to cations; (2) as potassium and 
sodium ions diffuse to equilibrium across the membrane, water enters 
the cell under the influence of the colloid osmotic pressure of the hemo- 
globin. This penetration of water results in swelling of the cell to a 
critical volume at which hemoglobin is released to the environment. The 
phenomenon of colloid osmotic pressure has been suggested to account 
for hemolysis following both ultraviolet (Leu et al, "42; Cook, 56) and 
X-radiation (Ting and Zirkle, 40) in mammalian red cells. The inhibi- 
tion of ultraviolet-induced hemolysis of frog red cells by sucrose 
(Maroney, ’59) lends support to the colloid osmotic hemolysis theory. 


For the colloid osmotic pressure of hemoglobin to be effective in 
causing hemolysis, the ions must be free to diffuse to equilibrium across 
the cell membrane. Energy is used to prevent the potassium and sodium 
ions from reaching diffusion equilibrium in the intact cell (Harris, *41; 
Danowski, 41; Maizels, ’54). Inhibition of metabolism is, then, one 
means for allowing ions to reach diffusion equilibrium. An example of 
this may be found in mammalian red cells kept at low temperatures 
(5°C). Potassium is lost from these cells and regained upon returning 
the cells to physiological temperature, provided that glucose is present 
(Harris, ’41). The results of the present study show that  ultraviolet- 
irradiated frog erythrocytes respire at the same rate as non-irradiated 
erythrocytes; yet cells exposed to ultraviolet radiation exhibited hemo- 
lysis while non-irradiated cells did not, It seems apparent that the cell's 
aerobic respiration is intact and an explanation other than the inhibi- 
tion of energy metabolism must be sought to explain hemolysis in the 
frog red cell. 
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The concept that metabolic energy is solely responsible for main- 
taining cation gradients has been criticized by Parpart and Green ('51). 
Rabbit red cells, when treated with n-butyl alcohol, gain sodium and 
lose potassium at comparable rates. Further, the recovery of potassium 
and sodium exchange after removal of the alcohol was much more rapid 
than would be expected if metabolism had been interfered with. They 
point out that it is more likely that the butyl alcohol has an effect on 
the membrane alone, permitting a rapid exchange of potassium and 
sodium. 


It is suggested, then, that ultraviolet radiation, in the case of the 
frog red cell, altered the cell membrane allowing potassium and 
sodium to diffuse freely. The cell then swelled because of the excess 
intracellular osmotic pressure. In this regard, frog erythrocytes were 
acting as did the mammalian red cell following ultraviolet radiation 
(Green, 56). It should be noted that the respiratory rates of irradiated 
and non-irradiated cells were the same in spite of the fact that up to 40 
per cent of the irradiated cells had hemolysed. This would indicate that 
the respiratory apparatus of the hemolysed cell was intact and func- 
tioning normally although further study of this phenomenon is desirable. 


SUMMARY 


Frog erythrocytes washed in phosphate buffered Ringer’s solution 
were exposed to ultraviolet light for various lengths of time. Respiration 
and hemolysis of both irradiated and non-irradiated cells were measured 
for about six hours following exposure to ultraviolet radiation, 


It was found that this radiation increases hemolysis without affect- 
ing oxygen consumption of these cells. 


The similarity of the mechanism of hemolysis in frog erythrocytes 
to the proposed mechanism of hemolysis in mammalian erythrocytes is 
discussed. 
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BALLOON-WING, A NEW MUTATION 
IN BLATTELLA GERMANICA (L.) 


Donald G. Cochran and Mary H. Ross 
Virginia Agriculture Experiment Station 


INTRODUCTION 


The study of spontaneous and induced mutations in such insects 
as Drosophila, Habrobracon, and several others, has been of great signi- 
ficance in the development of formal genetics. In fact, were it not for 
Drosophila, genetics might have progressed at a much slower rate. With 
regard to inheritance in the myriad of other insect species, however, 
there is still very little information available. For example, the mutation 
reported in the present paper is, to the writers’ knowledge, the first 
genetically controlled modification of the wing to be described in the 
ubiquitous German cockroach, Blattella germanica (L.). Although it 
might be argued that this is an extremely stable species and that muta- 
tions would, therefore, be uncommon, the facts do not seem to bear this 
out. During the course of the present investigation several other spon- 
taneously occurring variants have been found in the same inbred stock, 
Some of these traits are now under study, and the results will be 
published elsewhere. 


In Drosophila it is known that approximately 150 genes influence 
or modify the development of the wing (Emerson, 1955: 648). Muta- 
tions at several of these loci show some similarity to the one described 
here. In their compilation of known mutations in D. melanogaster, 
Bridges and Brehme (1944: 20, 23, 29, 26, 152, 95) listed the following 
traits: balloon (ba), Bubble (Bb), blister (bs), bloated (blo), puff 
(puf), inflated (if) and vesiculated (vs). Mutations of this type are 
also known in other species of Drosophila (Gordon et al., 1939: 51; 
Metz et al., 1923: 16) and in Musca domestica (Sullivan and Hirovoshi, 
1960: 214; Milani, 1956: 225). These mutations are all characterized by 
swellings on the wing, usually filled with fluid. 


The name “balloon” was chosen for the present mutation because 
the character seems to resemble most closely balloon-wing as described 
in the cheese skipper, Piophila casei, (Painter, 1917: 306) and, to a 
lesser degree, ballon-wing in D. melanogaster. Balloon in P. casei and 
in D. melanogaster is recessive and is not sex-linked. As will be shown, 
this is also true in B. germanica. 
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DESCRIPTION OF CHARACTER 


Balloon-wing roaches are characterized by a single, bubble-like 
swelling on each forewing and, usually, in each hindwing. During the 
course of this study 991 balloon-wing specimens were examined. On the 
forewings the swellings are oval or pear-shaped, usually from 0.24-0.28 
cm. in length and 0,15-0.17 cm. in width. In 85% of the specimens ex- 
amined the anterior ends of the swellings occurred opposite or proximal 
to the terminus of the first anal vein (1 A), but not within the anal 
area. This usual positioning of the swellings is shown in PI. I, A and B. 
A few swellings occurred within the anal area, and these were frequently 
associated with a minute swelling in a terminal or subterminal position 
(Pl. I, D). In certain crosses there was a tendency for the swellings to 
be more distal in position than usual, as shown in P1. I, C. Four speci- 
mens occurred with a swelling on one forewing only. 


A number of other variants of the forewing occurred with some 
frequency in certain crosses. Some of these are: divergent wings, wings 
with branching anal veins, wings shriveled or unfolded, dwarfed wings, 
and wings with one or more minute blisters. A swelling, filled with 
hemolymph, was also present on one or both sides of the pronotum in 
some balloon-wing roaches, but this character does not appear to be 
confined to particular crosses, 


On the hindwings the swellings are elongate, usually from 0.20- 
0.26 cm. in length and 0.10-,12 cm. in width. Sometimes the entire wing, 
with the exception of the anal area, is inflated. Hindwings of 869 bal- 
loon-wing roaches were examined. Of these, 818 had swellings on one 
or both wings. In 83% of the specimens the swellings occurred within 
the proximal two-thirds of the wing, but not within the anal area, as 
shown in Pl. II, A and B, Swellings within the anal area (P1. II, B) or 
in a terminal or subterminal position (Pl. II, C) were relatively rare, and 
confined to a few crosses. There were 103 roaches with one normal- 
appearing hindwing (98c%and 5 ), and approximately half this number 
(51, all o» ) which had no swellings on either hindwing. 


Several variants of the hindwing occurred which may be of genetic 
interest. As mentioned in the preceeding paragraph, two types of bubble- 
like swellings other than the usual balloon-wing were found. In addition, 
there were a number of tiny, blister-like swellings, a broken or incom- 
plete cubitus (Pl. II, C) and several other aberrations of the venation. 


GENETICS 
The balloon-wing roaches used for crossing came from a colony 


which apparently was breeding true for this trait. All individuals from 
the colony which were examined showed the balloon wing. They were 
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derived from those few mutants originally observed in a normal-appear- 
ing inbred stock. Those roaches involved in the crosses had the balloon in 
the usual position (P1, I, A and B) with a few in Pl. I, C. The wild 


type (normal) roaches came from a different stock. The following series 
of crosses were made: 


Series A. - Twenty-one crosses involving balloon-wing males and fe- 
males. 


Series B. — Twelve crosses reciprocally mating balloon-wing specimens 
with normal roaches. 


Series C. — Twenty-four crosses reciprocally mating balloon-wing speci- 
mens and normal-appearing hybrids. 


Series D. — Nineteen crosses involving hybrid males and females ob- 
tained from Series B. 


Series E. — Nine crosses reciprocally mating normal-appearing hybrids 
and roaches from the normal stock. 


Result of Crosses 


The results obtained from the five series of crosses listed above are 
presented in table 1. The reciprocal matings showed no significant dif- 
ference and are not listed separately in the table. The results indicate 
that the mutation under study is a simple Mendelian recessive. As would 
be expected, balloon x balloon (ser. A) produced all balloon offspring; 
balloon x normal (ser. B) produced all normal-appearing offspring; 
balloon x hybrid (ser. C) gave results which closely approximate a 
ratio of one balloon to one normal; hybrid x hybrid (ser. D) produced 
offspring in a ratio of three normal to one balloon; and hybrid x normal 
(ser. E) gave all normal appearing offspring. The results from series C 
and D were analyzed statistically by the chi square goodness of fit test. 
It was found that the chi square values were very low. This indicates 
that the deviations from a perfect 3:1 or 1:1 ratio are minor and are well 
within the range of deviations expected by chance alone. 


The four specimens mentioned in the preceding section which had 
swellings on only one forewing are believed to be the result of reduced 
expression of the trait. In order to test this hypothesis two of the speci- 
mens were successfully crossed, one with a normal and the other with a 
balloon-wing roach. In the F, of the former, all offspring were of normal 
appearance. In the F, there were 190 normal-appearing roaches and 
65 roaches with swellings on both forewings. This is also well within 
the expected chance deviations from a 3:1 ratio. In the second cross, 
there were 26 males and 32 females, all with the usual balloon-wing 
swelling on both forewings. Thus it would appear that in these two indi- 
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viduals, at least, the absence of a swelling on one forewing was a 
phenotypic variation. 


In B. germanica, as in other Orthoptera, the sex determination 
mechanism is of the XO type. If it is assumed that approximately half 
of the sperm carry the X chromosome, the sex ratio should be 1:1. In 
the present study an opportunity for checking this ratio was present. 
The results gave a total of 2132 males and 2085 females from 85 mat- 
ings. This is well within the expected range of variation for a 1:1 ratio, 


DISCUSSION 


In the present paper the character “balloon-wing” has been describ- 
ed and its inheritance elucidated. The associated problem of how the 
gene responsible for this trait causes the swelling to develop has not been 
investigated. However, the writers have observed the development of 
the swelling in several newly emerged adult roaches and found that it 
first shows as one, two or three tiny blisters which appear suddenly 
when the wing is partially unfolded. The blister or blisters enlarge 
quickly to form a swelling of the usual size by the time unfolding is 
complete. Normal-appearing veins are present before the appearance of 
the blisters and after formation of the swelling are still discernable on 
the lower, if not the upper, surface of the swelling. The only explana- 
tion for the trait which can be offered at this time is the speculative one 


that perhaps it is an inheritable weakness of the middle membrane of 
the wing. 


Along these same lines, it is interesting to note that the character 
(balloon-wing) does not always express itself fully in roaches. This is 
probably the result of variable expressivity as defined by Srb and Owen 
(1953:303). Evidence on this point comes from the two crosses reported 
in the previous section where the roaches showed a swelling on only one 
forewing. In both cases they were genetically of the balloon type. If it is 
assumed that this is the case, then a calculation of the per cent expression 
can be made, For this purpose it will be considered that complete ex- 
pression prevails in all those specimens in which a balloon appears on 
both wings of a pair. Accordingly, full expression of the trait in fore- 
wings is greater than 99 per cent, while in hindwings it is approximately 
80 per cent. This difference in expressivity in fore- and hind-wings has 
not been studied, but it may be genetically controlled. On the other 
hand, there is an explanation which does not require any genetic involve- 
ment. The hindwings normally lie in a folded condition between the 
dorsal body wall and the forewings. It is possible that the folds and the 
confined position confer enough structural reinforcement to the hind- 
wings that the swellings fail to appear in a higher percentage of hind- 
wings than forewings. 
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As adults, balloon-wing Drosophila lose the characteristic swelling 
and are best recognized by divergent wings and an extra plexus of veins 
(Waddington, 1940: 128). Balloon-wing roaches seldom have divergent 
wings and show no extra plexus of veins. The swelling on the wings 
are the only constant phenotypic expression of the mutation. 


SUMMARY 


A new mutation has been described for the German cockroach, 
Blattella germanica (L.). Based on analogy with similar traits in other 
insects, it has been named “balloon-wing.” The character has its pheno- 
typic expression as a swelling on each of the four wings of this insect. 
The swellings are relatively large, and on the forewings usually occur 
within the proximal half or near the center of the wing. On the hind- 
wings the usual expression is within the proximal two-thirds of the wing. 
In a small percentage of cases the swellings are unusually located. A 
number of other unusual variants, some of which are now under study, 
have also been described. 


In order to determine the mechanism by which this trait is inheri- 
ted, a series of crosses involving the balloon-wing stock and a normal 
stock have been made. The results showed that the swellings are the 
phenotypic expression of a simple Mendelian recessive gene. Four speci- 
mens (<1% of total) were found with a swelling on only one fore- 
wing. A number of specimens (approx. 20%), usually males, were 
found which lacked swellings on one or both hindwings. This is probably 
the result of variable expressivity, and, if so, the per cent of full expres- 
sion of the trait is considerably lower in the hindwings than in the fore- 
wings. A possible explanation for this difference together with a dis- 
cussion of several other points is included, 
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Table 1. — Results from a series of crosses involving a_balloon-wing 
and a normal stock of B. germanica. 


No. of normal-appearing No. of balloon-wing 


offspring offspring 

Series ¢ total o total 
A 

(balloon x balloon) 0 0 0 456 424 880 
B 

(balloon x normal) 417 423 840 0 0 0 
C 

(balloon x hybrid) 255 247 502 279 236 515 
D 

(hybrid x hybrid) 433 440 873 140 140 280 
E 

(hybrid x normal) 152 175 327 0 0 0 
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Explanation of Plate I 


Plate I. The drawings represent left forewings of B. germanica, and 
show the various expressions of the “balloon-wing” mutation, A and B, 
balloon in usual location on wing; C, balloon more distal than usual; 
D a small terminal swelling plus a balloon in the anal area. 


= 
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Explanation of Plate II 


Plate II. Similar to Plate I. Drawings of left hindwings. A, balloon in 
usual position; B. balloon in the usual position plus a balloon in the 


anal area; C, balloon more distal than usual plus an interrupted cubitus 
vein. 
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THE NEGATIVE BINOMIAL DISTRIBUTION: 
A REVIEW OF PROPERTIES AND APPLICATIONS* 


Joun J. BARTKO 
Virginia Polytechnic Institute 


1. JNTRODUCTION 


The negative binomial distribution (NBD) was discussed by Pascal 
and Fermat (Todhunter, 1865). However, the earliest known derivation 
and publication is due to Montmort in 1714. A particular form of the 
negative binomial, the inverse binomial sampling form, is sometimes 
referred to as the Pascal distribution. 


Many forms of the NBD appear throughout the literature. The most 
frequently occurring forms, as given in Table I, are discussed in this 


paper. 


TABLE I 
Forms of the NBD Parameters Comment 
(1) (** 1) r, p Pascal form; p+q=1, 
r is an interger. 
(2) r, p Pascal form; n=r-+x. 
k, p q=1+p; k not neces- 
(3) sarily an integer. 


The NBD, a two parameter distribution, has the same sample 
space (0, 1, 2, . . .) as the Poisson. In particular, the Pascal form (a 
term reserved specifically for the case where r is an integer) is a 
distribution over the non-negative integers, 


In more recent years NBD applications have increased, particularly 
in health and accident statistics, binomial sampling, population counts, 
data of contagious distributions and communications. 


2. ANALYTIC PROPERTIES 
Although many forms of the NBD are utilized, the Pascal form is of 
major importance. Its primary analytic properties follow. 


*Research supported by the National Institutes of Health and the Office of Ordnance Re- 
search, U.S. Army. 
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The characteristic function of the Pascal form is 

b(t) = p™[(1 - Accordingly, the mean is 

rqp7) and the variance rqp7*. In this distribution the variance 

is larger than the mean. From the form of p(t) it is easy to see 

that the sum of two independent negative binomial variates x, + X 
with the same p is also a negative binomial variate. 


The cumulant generating function is 
(t) = fn 


= -r - - 1)/pJ 


From the definition i.e., 


the factorial moments are given by = r(rtl).. »(rtk-1) 
A verv interesting and convenient property of the Pascal form is that 
its c.d.f. 
(* 
=O 


x I(r, yt1) 


(4) Fly) 


x 


can be computed by entering tables of the incomplete Beta function as 
was pointed out by Pearson and Fieller (1933). A presentation of their 
result follows. 


The remainder after y terms of the Maclaurin series in the integral 
form for a function g(x) is 


x 
R= ror | lat 
0 


Writing the Pascal form (1) as p*L(1-q)7"] and taking 
= (1-q)"" we have 


q 


Let (l-u) = (l-q)/(l-qtt)3 (5) with the factor p’ included 


then becomes 


9 
r 

Trt - 
iy y Ly Tau(1-q) 


fe) 
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and finally R = Ilys r), 


where I,(a, b) is the 
standard notation for the incomplete Beta function 


/pla,b). 


Therefore, 
and 
(6) F(y) =1 - r) = I(r, y+l) . 


It is interesting to note that the use of incomplete Beta tables for 
the evaluation of this c.d.f. has been rediscovered by Patil (1960). 


3. Limitinc Forms 
Three limiting forms of the NBD are discussed in this section, 
3.1 The Poisson Distribution. 
Assume the Pascal form (1) whose moment generating function 


(m.g.f.) is p’/(1-qe’)?. Let »% be fixed such that q = A/r, 
and lett then 


8 
(7) lim (—2-5)* er(e -1) 
l-qge 


the m.g.f. of the Poisson distribution, en ”*/x!, where x is the 
number of failures before the desired number of successes r. 


3.2 The Gamma Distribution. 
Consider the Pascal form (2) giving the distribution of sample size 


(ner) !(r-1)! ' 


n as 


Suppose each trial takes time At and that the probability of success p 
is proportional to dt, i-e., p= Adt. then time t = nit required 
to achieve r successes has the d.f. f(t) given by 
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= 1)1. tL 
(Se - 
f(t) = At)...(t - 
tr 


Let dt +0; then f(t) becomes the Gamma density 
ar e7 Mt 


8 f(t) = 
(8) I(r) 


which may be interpreted as the “holding time” before the rth occurrence 
of a rare event. 
3.3. The Logarithmic Distribution. 

In many applications of the NBD 


(9) PCa) (k-1) 


the number i observed in any sample may have all integral values 
including zero. However, in some cases the value i = 0 will not be 
observed, For example in the application of the NBD to the number of 
representatives of different species of butterflies obtained in a collection, 
only frequencies of numbers greater than zero will be observable, since 
by itself the collection gives no indication of the number of species 
which are not found in it. 


The following development of the logarithmic distribution is due to 
Fisher (1943). In nature the abundance of different species of the same 
group generally varies very greatly. Fisher found in a study of Malayan 
butterflies that the NBD which often fits such data well has a value of 
k so small as to be almost indistinguishable from zero. That it is not 
really zero for collections of wild species follows from the fact that the 
total number of species and therefore the total number not included 
in the collection is really finite. The real situation in which a large 
number of species are so rare that their chance of inclusion is small, is well 
represented by the limiting form of the NBD when k tends to zero. 


The limiting value k = OU cannot occur in cases where the frequency 
at zero is observable, for the distribution would then consist wholly of 
such cases. If, however, we put k = 0 in (9), write x for p/(1+p), so 
that x is a positive number less than unity varving with the size of the 
sample and replace the constant factor (k-1)! in the denominator by a 
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new constant factor a in the numerator, we have an expression for the 


expected number of species with i individuals, where i now cannot be 
zero, 


axt/ 
The total number of species expected is consequently 


(10) ax'/s = -a {n(1-x) 
i= 


the logarithmic distribution. The total number of individuals expected is 


@ 
(11) = ax/(l-a). 


Thus the two relationships (10) and (11) enable the series to be fitted 
to observational data, for if S is the number of species observed, and N 
the number of individuals the two equations 


(12) s=-afnll-x), N=ax/(l-x) 


are sufficient to determine the values of «@ and x. Tables which facilitate 


the solution of these equations and an example can also be found in 
Fisher (1943). 


4, Truncatep NBD 


In fitting data where the zero class is not observed it may be the 
case that a truncated NBD gives a better fit than the logarithmic 
distribution. Sampford (1955) considers the truncated NBD and gives 
procedures and tables to obtain iteratively either the moment or maximum 
likelihood estimates of the parameters. Included also is a discussion of the 
efficiency of these procedures. 


Simplified methods for the estimation of the parameters of the 
truncated distribution by a modified method of moments and by a 
modification of maximum likelihood procedures, due to Brass (1958) 


are discussed below. Brass also discusses the efficiency of these two 
procedures. 


Consider the NBD of the form (3) 


i 
13 p(i) = —(kti-1)'p 
(k-1) !i1(1+p)**2 


so that P(o) = 1/ (1+p)*. To obtain the corresponding proba- 
bilities (13) must be divided by [1 - P(o)}; writing w= 1/(1+p) 


d 

B 

c 

| 
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and it follows that 


i 
P(i) = (k+i-1)!y 
(k-1) 


where the subscript T denotes the truncated distribution. 


By considering the proportion in the first class of the truncated 
distribution and its mean and variance found by use of factorial moments 
Brass develops simple estimates for w, k and also for M, the mean of the 
complete NBD. The value M enters in his discussion of efficiency. Aside 
from M, the estimates are 


W 


=m(l - n,/n)/s* 
(14) _ Om - n,/n 
, 
where the bars denote the estimates. n. 


is the number of sample 
observations with measurement i and n is the total sample size. m and 


2 


s are the unbiased estimates of » and o, that is 
2 
con, (i m) 


Considering the log likelihood 


i 
L= £n.fn(% 


Brass obtains the following modified maximum likelihood expressions. 


k + n,/n 
(15) 
+ m 
m(k + n,/n) k + n,/n R_ R 
(16) = fn (— +4 + En, = 
k(m - n,/n) k +m j=l i=j 


where R is the largest observed value of i. The estimate from (14) gives 
a convenient first entry for K in the iterative procedure (16). w then 
follows immediately from (15). 

The above methods of fitting are illustrated on the data below also 
from Brass. The observations are the number of children ever born to a 
sample of mothers over forty years of age. 
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No. of Children Per Mother 
i 45 8&8 3 ii 


No. of 


m = 3.9912, s* = 5.9734, n,/n = 49/340 = 0.1441. 
Then from (14) ® = 0.572 and k = 5.00. 


For the maximum likelihood solutions write (16) as 


n,/n m-n,/n R R 
m=n,/n k k+m n j=l i=j i 
Then / / 
mn, /n m-n,/n R -2 R 


The iterative procedure is to calculate $ and $* for some value of k; 
a second value of k which will make  — closer to zero is then found by 
the aid of $+ on the assumption that the relation between $ and k is 
roughly linear over a short range. Further improved approximations are 
obtained by linear interpolation from the values of for the preceding 
two estimates of k. 


k or 
5.00 —0.00064 —0.00258 
4.75 +0.00008 

4.78 0.00000 


The first approximation to k is taken as 5.00 from the estimate by (14) 
The second approximation to k is 


—0.00064 
5.00 “—0.00258 4.75. 


One linear interpolation between 5.00 and 4.75 gives 4.78. From 
k = 4,78, @ can be found directly from (15) giving G@ = 0.561. The 
estimates again from (14) are k = 5.00 and © = 0.572. 


The estimates by the two procedures differ little. Brass gives the 
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efficiencies for both methods at about 95%. These efficiency values are 
read from tables which appear in his paper. The tables are based on 
estimates of k and M. In general the modified moments method of 


estimation is more valuable in that slight labor is involved and efficiency is 
high in many cases met in practice. 


5. OccurRRENCE 


Inverse binomial sampling is of primary importance in the considera- 
tion of the NBD. It occurs in the following general tvpe of situation. If 
a proportion p of individuals in a population possesses a certain 
characteristic, the number of observations in excess of r which must be 


taken to obtain exactly r individuals with the characteristic has the Pascal 
form of the NBD. 


Greenwood and Yule (1920) study the NBD in connection with a 
Poisson distribution with varying mean. 


Yule (1924), Furry (1939) and Kendall (1949) consider the dis- 
tribution in connection with birth-death processes. Yule is concerned with 
the mathematical theory of evolution; Furry with cosmic radiation and 
Kendall with the process in general. 


McKendrick (1914) and Quenouille (1949) discuss contagious dis- 
tributions. Irwin (1941) and Luders (1934) are concerned with the 
NBD as it occurs in multiple accident theory. 


The first of these situations, inverse binomial sampling, is the 
simplest mathematically and is the only one where the mathematical model 
is likely to hold exactly in practice. 


6. COMPUTATIONAL PROCEDURES ON THE PASCAL FORM 
6.1 Point Estimate of p. 
Occasions may arise in which the frequency of an attribute p is 
unknown. For the Pascal form the estimation of p from the sample data 
is a simple matter. 


Haldane (1943) considers the NBD (2) and derives an unbiased 
estimate of p, viz. 


(17) r-l 
n-l 


where n is the total number of items sampled in order to recover the 
predetermined number r. Haldane also develops an expression for the 
variance of this estimate. It is, however, approximate. Finney (1949) 
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derives an exact expression for the variance and also obtains an unbiased 
estimate of the variance. 


6.2 Interval Estimate of p. 


It may not have been generally realized that methods and tables for 
determining exact confidence limits on p in direct or binomial sampling 
may be adapted very easily to inverse binomial sampling. The proof 
which follows is due to Finney (1949). 


If py is the upper 4a confidence limit for a value of p estimated from 
inverse sampling, by definition it is the value of p for which the proba- 
bility of obtaining the quota r in n or more trials is exactly 4a. This 
value py is the same value of p which would give r-1 or Jess with the 
attribute in a sample of predetermined size n-1 in direct sampling. Like- 
wise p:, the corresponding lower limit, is the value of p which would 
give the quota in n trials or Iess with a probability of 4a, and is, there- 
tore, the same as the value of p which in direct sampling would give r or 
more with the attribute in a sample of size n. Thus, the rules for adapt- 
ing Table 41 of Bometrika Tables for Statisticians which gives confidence 
limits for the direct binomial can be stated: 


(i) For the upper confidence limit on p enter the table with 
= r-l, n = wl. 


(ii) For the lower confidence limit enter the table with c = 1, 
n=n. 


As an example of the use of these rules consider the following 
situation. An experimenter has concluded an inverse sampling plan in 
which he obtained four desired attributes in ten trials. Here r = 4, 
n = 10. To find the upper confidence limit on p enter the table with 
c = 3,n = 9. For the lower confidence limit enter with c —= 4, n = 10. 
Consequently, the two sided 95% confidence interval for p is 
0.11 < p < 0.70. The point estimate of p from (17) is 1/3. 


6.3 An Example of Inverse Binomial Sampling. 


Chapman (1952) considers an inverse binomial sampling plan in 
which it is desired to estimate the size N of a population, In Chapman's 
plan a number t of individuals are taken from this population. They are 
tagged and returned. In the second sample, n individuals are sampled in 
order to recover a predetermined number r of the tagged ones. 


If sampling is with replacement P(n) the probability of having to 
sample n individuals to obtain r tagged ones is given by 
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where r is the fixed parameter and n the random variable. €&(n) = rN/t 
and nt/r is an unbiased estimate of N. Chapman considers a normal 


approximation and a 2 approximation for confidence limits on N. 


However, a more expedient method is to use the rules stated in (6.2) and 
Biometrika Table 41. Let t/N = p, replacing N by nt/r, p becomes r/n. 


Then py and px can be determined such that Pip, <p Py? > l-a. 


On replacing p by t/N the confidence statement 
N t/P,) l-a is obtained. 


7. EsTIMATION OF PARAMETERS 


The estimation of parameters of other forms of the NBD is not as 
straightforward as it is in the case of the Pascal form. Parameter estima- 
tion in these cases is usually accomplished by the method of moments. 
However, when this method is not satisfactorily efficient, the method 
of maximum likelihood is commonly used. 


Consider first the method of moments. Fisher (1941) discusses the 


NBD (3) 


f(x) = 
(k-1) !x!{1+p) 


and finds the estimates of p and k by the method of the moments. These 


estimates are 
p* = - x)/x 


(18) and 


k* = x#/(s? - x) . 


Fisher also compares the determinants of the covariance matrices for 
large-sample estimates found by moments and by maximum likelihood 
obtaining an expression for the reciprocal of the efficiency. This further 
vields an expression, viz. 


(19) (1 + 1/p) (k + 2) 


which indicates when moment estimates are satisfactorily efficient, i.e., 
> 90%. This is the case when (19) exceeds 20 using the values of p 
and k found from (18). If the guide (19) indicates that moment 
estimates are not satisfactory then another method of estimation must be 
employed. 


Estimation by maximum likelihood is efficient. Haldane (1941) 
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considers the estimates of p and k by maximum likelihood. The estimate 
of p is 

(20) p = x/k 


where k is obtained iteratively from the following expression 


_ =n, Np 


+...+n 


n, is the observed frequency of x; , R is the the maximum value of x, 


R 
such that = Oh the total number 


R 
of observations and xX = = 


As an illustration of the two methods of estimation, consider the 
following data from Fisher (1941). The table gives a sample of sheep 
classified according to the number of ticks found on each. 


— < 123 4 5 6 7 8 9 10 
cay 7 9 8 13 8 5 4 3 0 1 2IN=60 
x = 3.25 = 5.9195, 


From (18) k = 3.75 and p = 0.821. With these values (19) is 
13.2 so that maximum likelihood estimates are the desirable ones. By 
(21) k = 3.75. p = 0.867 follows directly from (20). 


8. Hypotuesis TESTING 


Barnard’s (1946) method for testing the hypothesis Hj? Py = Pp in 
inverse binomial sampling is presented below. 

Consider the situation in which an instrument has been modified in the 
hope of removing a cause of failure. The sampling procedure consists in 
arranging to try the modified one, under similar conditions, and to con- 
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tinue to try each until a predetermined number r, of failures of the 
modified instrument and a predetermined number r, of failues of the 
unmodified instrument are observed (inverse binomial sampling). Then 
if p, refers to the modified instrument and p, to the unmodified one, and 
n,, n, are the number of observations made in the two cases, the proba- 
bility of the pair (nm, n,) is 

n,-1 n -1 r 

If p:. = p: = p, the probability of obtaining some pair of results having 
the same total n, + n, = N is 


Then the relative probability on the null hypothesis p, = p. of obtaining 
the pair (n,, n,) out of all results with he same total N is 


(N-(r,+r,) 1(n,-1) 1(n,-1) ! 


which is independent of p. Hence, if we sum over more extreme pairs 


having the same N, r, and r, we shall obtain the significance level of our 
result. 


For example, with r, = 1, r, = 2 if we find n, = 10, n.—4 we have 
the significance level 
2:11! 131! 1012! 
+ Tort * = VB 
which is non-significant at the 5% level. 


Barnard also develops an approximation procedure when p, and p, 
are small and n, and n, are large. The rule is to calculate 
R= (ny -r)/2)r/(n3-r,/2)r, and to enter the F tables with this value of 
R for (2r,, 2r,) degrees of freedom. 


For the case of a 2xf trial, provided p,, p.. . ., Pg are all 


small, Bartlett’s or Hartley's test can be used, Barnard did not elaborate 


this point. It seems worth-while, therefore, to examine how this situation 
can be referred to Bartlett’s test. 


Note that - and 
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- . Identify k, = 0? and 

(ay = ie., with 2r, degrees of freedom. 
Bartlett's test, m= (Nt - [)fn st ~ E{(nt-1)fn s?] where N* = total 
number of observation, n} = number of observations in the ith sample, 


{ = number of samples, then becomes in our notation - 


(23) M= 2ir,fn st - 2ir, fn s? 
where sf = - r,/2)/r;, i= 1,2, ...,f 


and - - r,/2)/tr,; ‘ 


To perform the test refer M to a 1° 


If “+++ = perform Hartley's test, max s*/min s’, and 
refer to Biometrika Table 31. 


9. y2 -TEsts FOR CONTINGENCY TABLES OF NEGATIVE 


BINOMIAL TYPES 


Steyn (1959) discusses the problem of the analysis of a contingency 
table where the frequency of one of the cells is predetermined. A 
summary of his work on this problem follows. 


Suppose sampling with replacement is performed in a_ population 
consisting of categories FE, FE, FE, and FE  in_ proportions 
Poi? Pio and Poo respectively, Individuals are drawn at random 
until r FE’s are found. Obviously the sample size is a random variable. 


The following table shows the observed frequencies for the negative 
fourfold table. This table is obtained when two events E and F can occur 
or fail either alone or together and sampling is continued until the rth 
joint failure has occurred, thereby fixing a priori the frequency in one 
cell but not the total in the sample. 


Y 
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Frequencies 
E E 
F ol 01 
* fo, * 44, * * 


Steyn develops the %* expression for the test on this table. The 
expression with one degree of freedom is 


where (p'/q') = (-(r-xty) + /(rtxty)” - Gxy)J/2r , 


(B/q) C-(rtx-y) + /(rtxty)* - yxy 


and (3/qqt) 


1 + (p'/qt) + (p/q) + (B/q)(B'/a")- 


The extension to the case where there are h possibilities for the event 
Eand f possibilities for F is also discussed. 


SUMMARY 


This paper is a review of the major literature dealing with the 
negative binomial distribution (NBD). An account of the forms of the 
distribution is presented. This includes limiting and truncated forms. 
Parameter estimation is treated in a few specific cases. The problem of 
testing P) = Pz» also that of testing p, = p, = -.- = Pf along 
with a discussion on contingency tables of negative binomial types is 
presented in some detail. 


The inverse binomial sampling form is emphasized throughout the 
paper. In the discussion of this form several extremely useful relations are 
discussed, Among these is a discussion of the use of the Tables of the 
Incomplete Beta Function for the evaluation of the cumulative distribu- 
tion function (c.d.f) Also discussed is the method of adapting Table 41 
of Biometrika Tables for Satisticians, giving confidence limits on p in the 
positive binomial, to the inverse binomial sampling case. 
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The opening of our new research 
center this month will provide us 
with an opportunity for even 
more extensive basic research 
into the physical sciences. 

The new laboratories and in- 
strumentation will facilitate 
comprehensive studies into the 
basic nature of tobacco, smoke 
and smoke flavor as well as in- 
tensify our corporate research 
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into packaging and adhesive ms 
terials. 

We are looking forward to the 
achievements to come from this 
new laboratory complex. Wé 
sincerely believe our greater ell- 
phasis on fundamental research 
will be an important factor 0 
the continued growth of bot) 
our company and the industries 


we serve. 
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NEWS AND NOTES 


(Editor's Notes — News contributions should be sent to the person whose 


name appears at the end of the appropriate sections. ) 


ENGINEERING SECTION 


Stricken by heart attacks during the Fall, both Professor Fred 
Bull of the chemical engineering department at Virginia Polytechnic 
Institute and Mr. Tilton E. Shelburne of the Virginia Council for 
Highway Investigation and Research are reported well on the road to 
recovery and expect to resume their duties early in 1961. 


Dean Lawrence R. Quarles of the School of Engineering, University 
of Virginia, attended a conference on University Reactors at Gatlinburg, 
Tennessee, sponsored by the Oak Ridge Institute of Nuclear Studies 
and presented a paper “Construction of the University Reactor.” In 
October Dean Quarles attended a meeting of the ANS-ASEE Committee 
on Objective Criteria for Nuclear Engineering Education in Chicago and 
also a meeting of the Region IV Executive Committee of the AIEE in 
Charlotte, North Carolina. Professor Orville R. Harris of the electrical 
engineering department of the University of Virginia attended a sympo- 
sium in New York in November on Magnetism and Magnetic Materials. 
In December Dr. Harris attended the annual meeting of the American 
Rocket Society in Washington. 


Professor R. E. L. Gildea spent a portion of the summer of 1960 on 
active duty with the U. S. Public Health Service during which time he 
helped organize a regional meeting on Radiological Health to be spon- 
sored during the spring by the University of Virginia, the Virginia De- 
partment of Health, and the U. S. Public Health Service. Major D. M. 
Crim of the civil engineering department at Virginia Military Institute 
ison leave of absence during the first semester of the 1960-61 academic 
year and is studying for a graduate degree in civil engineering at the 
to the University of Virginia. 

m this Professors Otis L. Updike and Robert M. Hubbard of the chemical 
x. We engineering department, University of Virginia, attended the annual 
er el- meeting of the A. I. Ch. E. in Washington in December. Professor 
search Updike was co-chairman for two sessions of the symposium on “The 
tor iD Unsteady State as an Analysis Tool in Chemical Engineering.” Professor 
> both Hubbard is chairman of two subcommittees of the A.I.Ch.E. Chemical 
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Engineering Education Projects Committee. A paper by Professor Hubbard 
“Teaching with a Pneumatic Analog” was published in the August and 
September issues of the ISA Journal. 


Professor Bull has announced that Mr. R. H. Pusey has assumed a 
teaching position on the staff of his department at VPI taking the place 
of Mr. G. F. Meehan who recently resigned, 


Professor Nelson F. Murphy of the chemical engineering depart- 
ment at VPI is a coauthor with Professor F, W. Bull of a bulletin from 
the Virginia Engineering Experiment Station entitled “Decolorization of 
Waste from Caustic Wash of Chlorine-Bleached Kraft Pulp.” Professor 
Murphy attended a meeting of the Committee on Education of the 
American Electroplaters’ Society in Hartford, Connecticut in February. 
In November he spoke to the Blue Ridge Branch of the American 
Electroplaters’ Society in Roanoke on the subject “How Acid Inhibitors 
Work.” 


Professor Robert Truitt of Virginia Polytechnic Institute has an- 
nounced that the name of his Department of Aeronautical Engineering 
has been changed to the Department of Aerospace Engineering. 


Mr. Tilton E, Shelburne of the Virginia Council for Highway Investi- 
gation and Research attended a meeting in Chicago in February of 
the Data Analysis Panel on the AASHO Road Test. 


Dean Lawrence R. Quarles of the Engineering School of the Uni- 
versity of Virginia attended a meeting of the American Society for 
Engineering Education-American Nuclear Society Committee on Ob- 
jective Criteria for Nuclear Education in Dallas, Texas in January. Also 
during that month Dean Quarles attended a national meeting of the 
American Institute of Electrical Engineers in New York. In February he 
attended a Conference on Higher Adult Education sponsored by the 
Southern Regional Education Board. He also met with the Regional 
Advisory Council on Nuclear Energy in Atlanta, Georgia in February. 


Professors Otis L. Updike and Robert M. Hubbard of the chemical 
engineering department of the University of Virginia attended the 
annual meeting of the American Institute of Chemical Engineers in 
Washington in December. Dr. Updike was co-chairman of two sessions 
of a symposium on the subject “The Unsteady State as a Tool for Analy- 
sis in Chemical Engineering.” Professor Updike has recently received a 
subcontract from the Oak Ridge National Laboratory for research on 
the dynamics of solvent extraction cascades. , 
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AGRICULTURAL SCIENCE 


Drs. D. G. Cochran and J. M. Grayson, V. P. I. Entomologists, have 
received a grant of $3,000 for studies in cross resistance and control of 
German cockroaches resistant to organic phosphate insecticides. The 
research will be in cooperation with the National Pest Control Associa- 
tion. Both investigators were speakers at recent Regional Pest Control 
Operator's Conferences, Dr. Cochran spoke at Purdue University and 
Dr. Grayson at the University of Massachusetts. 


Dr. Grover C. Smart, Jr., has joined the staff of the Holland sub- 
station as Nematologist. He received his B. A. and M. A. from the Uni- 
versity of Virginia and the Ph. D, degree from the University of Wis- 
consin. 


Several staff members in the Agronomy Department at Virginia 
Polytechnic Institute have recently received grant-in-aids. Dr. G. W. 
Thomas will study “Electrolyte Inhibition by Soil” on a $19,000 NSF 
grant, and Drs. R. E. Blaser and R. E. Schmidt have received $2,000 
from the American Potash Institute to conduct a project on the “Effect 
of Potash with Rates of Nitrogen and Phosphorus on Tifgreen Bermuda 
Grass Turf and Winter Growing Ryegrass.” Dr. Blaser has also received 
a $500 grant from the Golf Course Superintendents Association of 
American for an undergraduate scholarship in turf. 


New members of the Agronomy Department, Virginia Polytechnic 
Institute are Mr. H. H. Thomas, Mr. S. J. Ross, Jr., and Dr. J. Crews. 
Mr, Thomas and Mr. Ross will be stationed with the soil survey group 
in Madison and Dr. Crews will head up the department’s program in 
flue cured tobacco breeding. 


The Virginia Tech Chapter of A.D.S.A. (student Dairy Club) was 
honored as the outstanding dairy club in the United States at the recent 
A.D.S.A. meeting in Logan, Utah. 
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